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Abstract 

We present a classical integrable model of SU(N) isospin defined on complex 

projective phase space in the external magnetic field and solve it exactly 

by constructing the action-angle variables for the system. We quantize the 

system using the coherent state path integral method and obtain an exact 

expression for quantum mechanical propagator by solving the time-dependent 

Schrodinger equations. 
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There exist a number of integrable systems [l],[J which can be solved exactly. Classi- 
cal solvability is a question of finding the action-angle variables of the system. Quantum 
mechanically, one way of solving the system would be to find an exact expression for the 
propagator by solving the time-dependent Schrodinger equation. 

The classical models of point particles carrying isospin charge in interaction with external 
gauge field have been investigated for over two decades |§. Recently, attetion has been 
paid to them in relation to the non-Abelian Chern-Simons quantum mechanics [[§-§] which 
may have some application in condensed matter physics ||. In this Letter, we construct a 
completely integrable classical model of non-relativistic SU(N + l)(iV > 1) isospin particle 
interacting with constant external magnetic field and show that it can be solved exactly in 
the above way. We will be concerned with only isospin degrees of freedom and drop the 
spatial dependence all togather. 

Classical SU(N + 1) symmetry can be well described on complex projective space M. = 
CP(N) which is a symplectic manifold and therefore could be considered to be the phase 
space of classical mechanics. The symplectic structure is given by the Fubini-Study metric 
|TT| with a constant J 



where £ = (£i,---,6v) is the coordinates of CP(N) and |£| 2 = E^ =1 C m Cm = We 
introduce the notation £ A = £ n ) (m, n = 1, • • • , N) and write the symplectic two form 
Eq.(||D as tt = ~tt AB d£ A Ad£ B . Then, the Poisson bracket between F, H E C°°{M) is defined 



VI = 2iJ 



i + iei 2 (i + iei 2 ) 2 



(i) 



by 



{F, H} = Q AB d A Fd B H 



(2) 



where £1 



is the inverse matrix of VLab- The use of Eq.(|l]) results in the following expression 



for Poisson bracket: 




(3) 
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where g mn is the inverse of the Fubini-Study metric given by 

9 mn = ^ja + m(s mn +c m u. (4) 

Now, the group SU(N+1) acts transitively on CP(N) and leaves the symplectic two form 
Eq.(|l]) invariant. It can be shown || that the group action g = exp(itT a ) e SU(N+l),t G R 
generates the following Hamiltonian vector field on M. 

X a = i [(T a ) mo + (T a ) mn £„ - (T a ) 00 U - {T a )UnU J- + (cc). (5) 

Here T a 's form the generators of the Lie algebra Q = su(N + 1): 

[T a ,T b ] =if abc T c (6) 

We use the normalization Tr(T a T b ) = {l/2)5 ab and the notation that T a is (iV + 1) x (iV + 1) 
matrix with element {T a )ij = 0, 1, • ■ ■ , N). The above Hamiltonian vector field defines 
the isospin function (so-called momentum mapping function 0]) Q a on A4 by 

X a \tt + dQ a = Q. (7) 

One can show that Q a can be expressed as follows |],|| 

N 

Q a = - 2J E <(n^,-L o _ , _ uoC • (8) 

i,3=0 Vi+ICI 2 

It can be easily checked that X^s form su(N + 1) Lie algebra [X , X&] = f abc X c and Q a 
satisfy su(N + 1) Poisson-Lie algebra : 

{Q a , Q b } = f abc Q c - (9) 

To define our integrable model with the use of isospin functions Eq.(|8|) and @, we need 
an explicit representation for the generator T a 's. Recall that T a 's can be decomposed into a 
maximally commuting set {H m }(m = 1, ■ • • , N) which is called the Cartan subalgebra and 
a set of ladder operators E a one for each root vectors a. We use the defining representation 



of SU(N + 1) where the diagonal H m 's are defined as follows [p|: H m has m l's along 



the diagonal from the upper left-hand corner. The next diagonal element is — m to make 
it traceless. The rest of the diagonal element (if any) are zero. So with the normalization 
defined previously, we have 



m— 1 



J2 dik^jk ~ m5i, m 5j,m)/ J2m(m + 1) 



(10) 



fc=0 



The vector field generated by exp(itH m ) can be expressed using Eq.(|5|) as 



A. 



—% 



d N d 



+ (c.c) 



(11) 



'2m(m + 1) 

In the above, the index m is not summed. The corresponding isospin functions Eq.(^) are 
given by 

-23 



m—l 



Q' 



,Um=U0t,m 



(12) 



'2m(m + 1) fc=o 

We define an integrable model on M. with the Hamiltonian given by a linear combination 
of the above isospin functions 



/V 



(13) 



m=l 



with real constants fi m . When restricted to SU(2) case with N = 1, this model describes a 
classical (iso)spin in the constant external magnetic field as can be readily checked in terms 
of stereo graphical projection £ = tan(#/2)e~^: 



H = fiQ, Q = —3 cos 9. 



(14) 



with the magnetic field given by fi and 3 represents the magnitude of classical (iso)spin. 
By construction, there are N conserved quantities Q m which are in involution 



{H,Q m } = 0, {Q m ,Q n } = 



(15) 



and the system is completely integrable |]|§. The Liouville's theorem states that the man- 
ifold M.q defined by the level set of the functions Q m , 



M Q = {x:Q m (x)=q m ,x = (Z,0} 



(16) 



with g m =constant, is a smooth iV-dimensional manifold diffeomorphic to the torus 

T N = {(<f )u --- j( f )N )mod 2tt}. (17) 

and invariant with respect to the phase flow with the Hamiltonian. In this case, we can 
find the action variables (Ji, • • • , Jjv) conjugate to the angle (0i, • • • , 4>n) so that the original 
symplectic structure Eq.(Q) is expressed by the canonical two form 

N 

fi=^d/ m Arf0 m . (18) 

m=l 

The explicit form of the action variables can be found by the use of stereo graphical 
projection. Let us introduce the polar angles (9i, • • ■ , 6*tv)(0 < 6 < ir) and (0i, • • • , 0at) via 



£i = tan(0i/2) cos(0 2 /2)e~ 
£ 2 = tan(0i/2) sin(0 2 /2) cos(0 3 /2; 



6v-i = tan(0i/2) sin(0 2 /2) • • • sin(0jv_i/2) cos{8 N /2)e 
6v = tan(0i/2) sin(0 2 /2) • • • sin(0 JV _i/2) mn(6 N /2)e 



-i</>N 



(19) 



It is interesting to find that the action variables I m are given by the following formula: 



I m = 2Jsin 2 (# 1 /2) • • -sin 2 (0 m /2) cos 2 (0 m+ i/2) (m < N) 
I N = 2Jsin 2 (#i/2)- ••sin 2 (^_ 1 /2)sin 2 (^/2) 



(20) 



It can be expressed more compactly as I m = 2J£^£ m /(l + |£| 2 ). One can check that substi- 
tution of Eq.fllD into Eq.(0) prod uces Eq. fll8|) . The above action variables enables Q m to 
be expressed as a linear combination of them: 



w2m(m + 1) 
Substituting into Eq.(O), we obtain 



N 

(m + l)J m + Yl h~2J 

fc=m+l 



(21) 



N 



H ^nJn + H 



(22) 



n=l 



where Ho is a constant and the Larmor frequency u n is given by 



N Urn 

UJ n = j m = K m + l ¥nm + Onm] (23) 

m=i J2m(m + 1) 

with the step symbol 9 nm , 

{0, n < m 
(24) 
1, n > m. 

With Eq. (p2|) , we have completely solved the system and the solution is given by 

m (t) = 9 m (Q), 4> m {t) = 0m(O) + 0J m t. (25) 

Instead of performing the path integral in terms of the above Darboux variables |12,13 



we quantize the above classical model using the coherent state path integral method 
|i4Hl6| and obtain the exact quantum mechanical propagator by solving the time-dependent 
Schrodinger equations. We first construct coherent states on CP(N). Let us consider |0), 
the highest weight state annihilated by all positive roots of SU(N + 1) algebra in Car- 
tan basis. Then for CP(N) with given P = 2J (P 6 Z + ) we have an irreducible repre- 
sentation (P, 0, • ■ ■ , 0) of SU (N + 1) group |F7| and there are precisely N negative roots 
E a , a = 1, 2, ■ • • , iV such that E a \0) ^ |0). Let us label {E a } = {E m }. We define a coherent 
state on CP(N) corresponding to the point £ = (£i, • • • , 6v) by [P||l6[| 

|P,e) = exp(^e mJ B m )|0) (26) 

m 

Notice that this definition differs from the usual one by the normalization factor. We have 
chosen this definition here because in the subsequent analysis, £ and £ can be treated in- 
dependently and the overspecification problem can be side-stepped |18|-p0[]. We denote 
\P, £) = |£) from now on. The coherent states which we have defined on CP(N) have the 
following two properties which are essential in the path integral formulation. One is the 
resolution of unity, 



Ml*) * \ILj =I i (27) 
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where = cd£d£/(l + \£\ 2 ) N+1 with a constant c is the Liouville measure [1161. The 

other is reproducing kernel, 



We are interested in evaluating the propagator 



G(f,e , ;t) = (rie- iffi ie / ) 



(28) 



Inserting Eq.fl27|) and using Eq.( pgf) repeatedly, we obtain the following expression 



(29) 



G(Z",?;t) = I D/iexp|2Jlog(l + ^(t/))+« jf dt 



2J& 



(30) 



i + iei 2 

where H(£,£) = * s the classical Hamiltonian given by the Eq.([l3|). The 

boundary conditions in the path integral is given by = £' and £(£/) = We introduced 
which is only a superfluous variable because the result of path integral Eq.(|3~0"D does 
not depend on this variable. It depends only on £" and 
The equations of motion 



(31) 



are those of collections of N harmonic oscillator although the Hamiltonian Eqs.fll3|) and ([12] 
appears to be highly nonlinear: 



Cm ~ i^mCm = 0, £ m + iw m £ m = 0. 



Here u; m is given by Eq . (|23|) . The solutions are given by 



(32) 



«/ iu m (t-t") 



' e —iu>m{t—t') 



(33) 



Using the above fact, we can evaluate the propagator by the semiclassical approximation 



method |ZT| and the result agrees J22] with the exact propagator which is obtained by solving 
the time-dependent Schrodinger equation as is adopted in this paper. 

In order to have explicit operator form for the Hamiltonian to set up the Schrodinger 



equation, we resort to geometric quantization | 23fl . In the geometric quantization of classical 



phase space M. = CP(N) with symplectic structure Q, we quantize classical observables 
F = F(Q a ) which are functions of only Q a 's satisfying the Poisson-Lie algebra Eq.(|9[). The 
prequantum operator corresponding Q a is given by 

Q a = -tX a -X a \Q + Q a (34) 

where is the canonical one form Q = dQ given by 

Since our Hamiltonian Eq. ([T3|) is linear in Q a 's, there is no normal ordering ambiguity. Also 
the polarization is chosen such that £?(£",£';£) is a function of £" but not of Hence we 
get the following time- dependent Schrodinger equation for the propagator: 



d - N 

m =i J2m(m + 1) 



dt 



(m + l)e. 



d 



N 



m ftp*!! 



EC*/' ^ 
>fc=m+l UC *k 



(36) 



We note that the above equation is invariant under the gauge transformation of the canonical 
one form by — > + dA(£, ^). The boundary condition is given by 



"C'\2J 



(37) 



It is remarkable that the solution can be obtained in a closed simple expression by 

/ N \ 2J N 

G(C", 1+E CC exp{-z^ m t} exp{z J ^ n m /y/m(m+ l)/2t}. (38) 

V m=l / m=l 

The above results reduces to the exact propagator and reproduces the Weyl character formula 
for 517(2) case . 



In conclusion, we presented a classical integrable model of SU (N) isospin in interaction 
with external constant magnetic field neglecting the spatial dependence completely and 
solved it exactly both classically and quantum mechanically Especially, we obtained the 
exact expression for the quantum mechanical propagator by solving the time- dependent 
Schrodinger equation set up by the method of geometric quantization. This method could 



S 



be applied to other integrable models such as the one in which the Hamiltonian is a sum of 
quadratic functions of Q m of Eq.(Jf2]) and other coadjoint orbits |I7[ of Lie group including 



the non-compact case. These and related issues will be reported elsewhere . 
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